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基于 规则 网 格 的 层次 化 布料 动态 模拟 方法 


郑 “ 华 ， 王 美丽 ， 李 书 琴 ; 


(西北 农林 科技 大 学 ,信息 工程 学院 , 陕西 咸阳 712100) 


摘 要 : 在 布料 动态 模拟 仿真 过 程 中 ， 收 化 速度 和 模拟 效率 是 两 个 核心 指标 ， 可 以 很 好 地 反映 布料 在 动态 过 程 中 的 模 
拟 效 果 。 针 对 布料 动态 模拟 中 的 收敛 速度 慢 、 模 拟 效率 低 的 问题 ， 提 出 了 一 种 基于 规则 网 格 的 层次 化 模拟 方法 ， 实 现 
oo ned aml 利用 层次 化 思想 在 原始 网 格 的 基础 上 构建 层次 化 约束 网 格 。 在 这 个 过 程 
中 ， 可 以 采用 不 同 的 决策 函数 对 网 格 进行 精简 ， 构 造 出 更 加 满足 目标 要 求 的 约束 网 格 ， 在 构造 完成 后 利用 提出 的 权 值 


关联 模型 对 各 层次 进行 再 矫正 。 在 模拟 过 程 中 ， 原 始 网 格 利 用 层次 化 约束 网 格 从 最 粗 层 到 最 精细 层 进行 收敛 矫正 ， 


有 


效 避 免 了 传统 为 提高 收 你 速度 而 单一 的 增加 约束 矫正 选 代 次 数 的 问题 。 为 了 检验 模拟 性 能 ， 布 料 在 周期 钟 摆 风 场 下 进 
行 了 实验 。 实 验 表 明 ， 在 基于 规则 网 格 的 层次 化 方法 模拟 下 ， 能 很 好 地 解决 传统 模拟 的 约束 震荡 问题 ， 并 且 效率 高 、 


稳定 性 好 。 
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Cloth simulation based on grid-hierarchical dynamic simulation method 


Zheng Hua, Wang Meili, Li Shuqin’ 
(College of Information Engineering, Northwest A&F University, Xianyang Shaanxi 712100, China) 


Abstract: In the process of cloth dynamic simulation, convergence and efficiency are two core indicators, which can well reflect 


the performance of simulation. To approve the convergence and quality in cloth dynamic simulation, this paper proposed an 


approach which based on grid-hierarchical simulation. In this method, this paper constructed a hierarchical constrained grid on 


the basis of the original grid using the hierarchical idea, in this process, various decision functions could be used to streamline 


the grid to adapt different target requirements. After the construction was completed, it used the proposed weighted association 


model to rectify each level. In the process of simulation, the proposed approach used hierarchical constraint grid to constrain the 


original grid from the coarse layer to the finest layer to improve the convergence speed, which effectively avoided the problem 


of increasing the number of iterations in the traditional way. To show the simulation performance, it simulates dynamic cloth 


under the cycle pendulum wind field. The results shows that under the simulation of grid-hierarchical method, it can well solve 


the constraint oscillation problem of the traditional simulation, and the proposed method is more efficiency and stable. 
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对 于 布料 模拟 的 建 模 有 以 下 三 种 比较 典型 的 方法 : 


0 引言 


第 一 种 


就 是 根据 一 些 经 验 性 的 几 人 


1 


方程 来 表示 材质 特有 的 物 开 
在 在 计算 机 动画 中 ,布料 动画 模拟 是 重要 的 研究 方向 之 一 ， ”其 优点 是 可 以 高 效 生成 的 良好 的 模拟 效果 。 典 型 的 模型 有 悬 链 


特性 ， 


受到 了 很 多 研究 者 的 关注 。 如 何 让 一 个 平面 真实 的 反映 出 其 自 ” ” 线 模型 、 纯 几 


, 洒 


变换 形变 模拟 模型 、 基 于 纹理 生成 皱 裙 模型 与 


身 的 连贯 性 、 拉 伸 性 、 磁 撞 性 和 受到 外 力 的 反应 成 为 了 研究 的 。 双 层 几何 模型 (衣服 层 、 皮 肤 层 )。 第 三 种 是 基于 布料 的 摩擦 、 
重点 I。 目 前 已 有 的 众多 算法 在 效果 、 效 率 以 及 真实 性 上 的 差 弹性 系数 等 物理 特征 ， 通 过 构造 三 角 网 格 或 矩形 网 格 的 结构 力 


异性 较 大 ， 所 以 研究 一 个 综合 性 更 强 、 效 果 更 好 、 效 率 更 高 的 。 ”学 模型 和 能 量 状 态 模 型 ， 然 后 根据 动力 学 方程 进行 数值 求解 ， 


模拟 算法 是 当前 布料 模拟 的 重点 研究 方向 ， 同 时 也 为 后 续 其 他 ”其 优点 是 可 以 更 加 逼真 的 显示 布料 。 该 方法 又 分 为 连续 体 模 型 
材质 的 模拟 葛 定 基础 。 布 料 动画 模拟 包括 了 对 于 布料 建 模 、 数 。 和 离散 体 模型 。 连 续 体 模型 典型 的 有 弹性 形变 模型 、 线 性 应 变 
值 计 算 与 模拟 方法 三 个 方面 。 力 模型 、 非 线性 应 变 力 模 型 和 细 分 有 限 元 模型 。 目 前 对 于 离散 
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型 模型 研究 最 为 深入 ， 典 型 的 代表 有 Breen 等 人 外 提出 的 粒子 。 束 网 格 的 质量 严重 影响 了 收敛 的 速度 和 模拟 效果 ， 并 且 在 多 层 
模型 、ProvotD] 提 出 的 质点 -弹簧 模 型 、Kikuuwe 中 提出 的 通用 ”约束 网 格 的 层次 关联 上 存在 父子 粒子 关系 不 确定 的 问题 。 在 X- 
粒子 模型 。 第 三 种 就 是 综合 以 上 两 种 方法 的 优点 ， 即 几何 模型  PBD09 方 法 中 , 虽然 提出 了 新 的 控制 布料 模拟 的 手段 , 但 是 对 于 
的 高 效 性 和 物理 模型 的 逼真 性 。 除 了 以 上 几 种 典型 的 模型 外 ， 层次 化 模拟 的 不 确定 问题 、 模 拟 时 存在 的 局 部 收敛 问题 并 没有 
Lu 等 人 中 提出 了 一 种 基于 NURBS 的 物理 模型 。 该 模型 在 布料 仿 。 解决 。 在 最 近 提 出 的 基于 采样 的 层次 化 模拟 方法 中 中， 虽然 能 


真 程度 上 达到 了 新 的 高 度 ， 但 是 仍然 存在 模拟 效率 低 和 布料 收 ” 较 好 地 解决 层次 化 的 网 格 质 量 问 题 ， 但 是 网 格 质量 的 不 可 控 性 
敛 速度 慢 等 问题 。 使 得 其 在 应 用 中 和 理论 上 存在 着 很 大 的 不 确定 。 
对 于 数值 计算 有 以 下 四 种 典型 的 求解 模型 :a) 显 式 Eluer 为 了 解决 以 上 问题 ， 本 文 提 出 了 一 种 新 的 基于 位 置 的 层次 


积分 方法 , 该 方法 相 较 于 其 他 方法 算法 效率 更 高 ; b) 隐 式 Eluer ”化 动态 算法 。 本 文 方法 通过 自 定义 决策 函数 来 构造 新 的 约束 网 
积分 方法 ， 主 要 应 用 于 解决 刚性 材料 的 模拟 ，c)4 阶 Runge- ，” 格 ， 将 约束 网 格 的 构建 与 决策 函数 挂钩 ， 不 再 单纯 地 依赖 具体 
Kutta 积分 框架 , 该 方法 的 优点 是 模拟 精度 高 、 模 拟 效果 良好 ， 算法 ， 在 将 非 线性 Gauss-Seidel 求解 器 转换 为 基于 非 线 性 多 
缺点 是 算法 效率 较 低 , 不 能 实时 显示 结果 ; d) Verlet 积分 框架 ， ” 重 网 格 的 过 程 中 ， 提 出 了 一 种 新 的 权重 关联 方式 ， 使 得 层次 间 
该 方法 针对 几 种 方法 的 缺点 进行 了 改进 ， 耗 时 少 且 模 拟 效 果 精 ”的 关联 不 再 依赖 于 不 确定 个 数 的 粒子 群 ， 丛 而 代 之 的 是 具有 三 
度 较 高 。 本 文 在 布料 建 模 上 使 用 通用 粒子 模型 ， 在 计算 方法 上 。” 角形 约束 的 父子 关系 群 ， 使 得 约束 矫正 效果 更 加 明显 ， 并 且 解 
采用 Verlet 方法 对 目标 网 格 进行 求解 。 决 了 迭代 和 时 间 步 长 相关 的 刚度 问题 。 
在 布料 模拟 的 过 程 中 ， 力 学 模型 中 传统 的 方法 是 利用 牛顿 为 了 展示 布料 动态 模拟 中 形变 较 大 的 位 移 和 角度 变化 的 
力学 方程 求解 ， 通 过 利用 加 速度 与 质点 质量 对 动态 物体 的 运动 效果， 展示 布料 在 动态 过 程 中 的 整体 稳定 性 ， 需 要 建立 基于 实 
轨迹 进行 求解 ， 其 中 存在 与 真实 相差 较 大 的 模拟 误差 。 因 此 ， ” 时 风 场 下 的 布料 模型 。 在 现实 场景 中 布料 的 形变 往往 伴随 着 风 
基于 此 方法 研究 的 主要 目的 是 尽 可 能 降低 误差 ， 以 此 提高 模拟 。” 场 的 变化 3， 周期 性 飘动 成 为 布料 的 典型 特征 中 。 为 了 更 好 地 
效果 ， 而 巨额 的 运算 导致 了 显示 模拟 结果 时 出 现 卡 顿 和 逼真 性 ”对 布料 进行 仿真 测试 ， 笔 者 将 风 场 环境 应 用 到 了 模拟 测试 。 
下 降 等 问题 。 、 

在 碰撞 检测 方面 上 ， 传 统 方法 需要 依赖 于 包围 体 结构 进行 1 理论 基础 
加 速 检 测 ， 如 常见 的 方法 如 球 包围 、 规 则 算 形 包围 、8-DOP 等 ， ”1.1 布料 建 模 
虽然 可 以 检测 碰撞 ， 但 是 明显 影响 了 模拟 的 效果 。 虽 然 近 些 年 本 文 使 用 目前 流行 的 弹簧 一 质点 模型 为 布料 建 模 ， 其 模型 
也 有 一 些 优秀 的 碰撞 检测 算法 如 TightCCDG、 虚 拟 分 解 方法 四 图 示 与 说 明 请 参见 文献 [13]。 对 于 模拟 动态 物体 , 首先 需要 包括 
等 使 用 GPU 获得 了 在 碰撞 检测 方面 上 的 提速 , 但 在 基础 模拟 算 ”模型 的 原始 位 置 、 各 粒子 的 动态 速度 与 各 粒子 的 质量 等 基本 物 
法 上 的 效率 依旧 是 仿真 效果 存在 的 一 个 瓶颈 。 而 基于 位 置 的 动 理 信息 。 其 次 需要 建立 各 粒子 间 的 约束 关系 。 其 中 本 文系 统 使 
态 方法 -上 能 较 好 地 解决 这 些 问题 。 
基于 位 置 的 动态 模拟 (PBD) 是 游戏 和 交互 式 应 用 中 变形 
体 的 实时 仿真 的 流行 方法 。 该 方法 的 简单 性 和 和 鲁 棒 性 特别 具有 E ee 
EE 0 0 As 
随 着 知名 度 的 提高 , PBD 的 局 限 性 也 越 来 越 被 关注 。 PBD 的 行为 。” 息 ， 即 各 粒子 间 的 拉 伸 范 围 以 及 剪 切 力 等 限制 ， 使 用 C 表示 。 
取决 于 模拟 的 时 间 步 长 和 大 代 次 数 中 。 具体 而 言 ， 随 着 迭代 次 。 1.2 受 力 分 析 
数 的 增加 ， 或 随 着 时 间 步 长 的 减少 ， 约 束 条 件 变 得 非常 死板 。 在 弹簧 质点 模型 中 ， 质 点 的 受 力 有 两 种 ， 一 种 是 由 质点 间 
当 创建 具有 各 种 材料 类 型 的 场景 时 ， 参 数 的 这 种 耦合 是 特别 可 ”的 弯曲 弹 徐 、 剪 切 弹 得 对 本 身 所 产生 的 内 力 ， 另 一 种 则 是 由 重 
能 引发 问题 的 , 如 与 近乎 刚体 相互 作用 的 柔 体 ,在 这 种 情况 下 ， ” 力 和 其 他 仿 自然 力 组 成 的 外 部 力 。 所 以 布料 的 受 力 可 如 式 (1) 
提高 迭代 次 数 来 获得 一 个 对 象 的 刚度 可 能 会 无 意 中 改 变 其 他 仿 。 ”所 示 。 
真 对 象 的 行为 。 这 往往 需要 在 整体 范围 内 重新 调整 刚度 系数 E 可 
使 得 创建 出 重复 使 用 的 仿真 资产 是 非常 困难 的 。 例 如 ， 在 布料 0 
模型 中 设置 拉 伸 和 弯曲 约束 的 相对 刚度 ， 就 可 以 单独 设置 迄 代 
次 数 。 更 糟糕 的 是 ， 选 代 次 数 的 影响 是 非 线 性 的 ， 难 以 直观 地 
调整 参数 ,或 者 很 难 将 重新 定义 的 值 作为 沈 代 计数 的 简单 函数 。 ” 示 质 点 所 受 的 外 力 。 对 于 受 力 的 具体 求解 请 参见 文献 [14]。 
1 基于 位 置 的 层次 化 动态 模拟 方法 将 (PBD) 的 非 线性 Gauss- 
Seidel 求解 器 转换 为 基于 非 线性 多 重 网 格 的 算法 , 从 而 显 着 增 
加 了 PBD 方法 的 收敛 速 度 ， 同 时 又 保持 了 各 自 的 特点 。 但 是 在 
具体 的 实现 中 笔者 发 现 ， 在 构造 多 重 网 格 〈 约 束 网 格 ) 时 ， 约 


f 诬 


二 


用 x 表示 布料 网 格 中 粒子 i 的 位 置信 息 ， 相 对 应 的 有 该 粒子 的 


其 中 ， 厂 ,表示 质点 所 受 的 力 ， 厂 , 表示 自身 的 内 力 ， 厂 , 则 表 


1.3 ”约束 矫正 
在 布料 动态 模拟 中 , 距离 约束 可 以 理解 为 质点 间 的 理论 弹 
得 长 度 ， 即 质点 1 与 了 间 的 理论 距离 。 如 果 质 点 间 的 实际 距离 
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不 等 于 理论 距离 ， 则 此 时 就 需要 对 质点 与 j 进行 一 定 的 位 移 
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于 原始 距离 。 当 函数 值 不 为 0 时 ， 就 需要 利用 约束 方程 进行 矫 


矫正 。 换 名 话说， 就 是 关于 距离 的 约束 矫正 。 不 同 的 约束 中 需 
要 不 同 的 约束 矫正 方法 ， 基 于 位 置 的 动态 模拟 则 提供 了 如 何 进 
行 约束 矫正 的 方法 。 
1.4 基于 位 置 的 动态 模拟 

基于 位 置 的 动态 模拟 (PBD) 09 是 一 种 可 以 直接 调控 三 维 
模型 的 位 置 坐标 的 方法 ， 它 以 友 代 Gauss-Seidel 方式 解决 位 
级 别 的 约束 问题 。 使 用 PBD 算法 模拟 可 以 更 方便 地 处 理 模拟 
过 程 中 所 需 添加 的 约束 ， 如 布料 端点 互相 的 距离 约束 、 布 料 的 
弯曲 约束 或 者 是 外 界 的 物体 碰撞 约束 等 。 

在 PBD 算法 中 通过 设 定 约束 函数 的 方法 为 模型 添加 常规 
约束 ， 而 不 是 仅仅 通过 计算 动能 传导 和 常规 力 的 方法 来 约束 模 
型 。 在 该 方法 中 ， 首 先 使 用 动力 学 方程 和 数值 计算 方法 提供 位 
变化 预测 ， 对 模型 所 有 点 的 坐标 根据 所 受 的 力 通过 动力 学 规 
律 得 到 下 一 个 时 间 周 期 的 位 置 坐标 ， 然 后 通过 模型 位 置 坐标 来 
检测 此 时 的 形状 、 位 置 是 否 满足 了 之 前 添加 的 常规 约束 ， 如 果 
不 满足 ， 则 根据 约束 矫正 方程 修改 点 的 坐标 ， 将 点 移动 到 满足 
约束 的 位 置 上 。 具 体 如 下 算法 1 所 示 。 

算法 1 基于 位 置 的 动态 模拟 管线 
输入 : 布料 网 格 ， 外 界 环境 参数 (如 重力 、 风 力 、 湿 度 )。 
输出 : 布料 动态 运动 轨迹 。 
1. for all 粒子 i do 
初始 化 x; 一 各, 二 Vw, = 二 1/m, 


i 


end for 
loop 
for all 粒子 i do 
Vv, vt+At:w,:f,,() 
end for 
for all 粒子 i do 
9. Pp; Xt+At:y, 
10. end for 
11. loop 从 代 次 数 
12. 约束 矫正 
13. end loop 
14. for all 粒子 i do 


15. Vv; €(p,—X)/At 


名 ”人 


16. Xp; 
17. end for 
18. end loop 


在 该 方法 中 通常 使 用 约束 方程 对 弹簧 质点 间 的 距离 进行 
矫正 ， 首 先 通过 牛顿 动力 学 方程 对 质点 的 位 置 进 行 预测 ， 然 后 
使 用 方程 对 超出 预期 赣 值 的 质点 进行 矫正 ， 使 其 恢复 到 理论 位 
内 ， 最 后 在 将 这 些 矫 正 后 的 质点 演 染 模拟 。 对 于 距离 约束 使 
用 的 方程 如 式 (2〉 所 示 。 


Cpi,p») pi -ps|-d (2) 


其 中 : C 中 的 pl 、ps 表示 要 操作 的 两 个 粒子 ，| pi 一 p, | 表 


示 两 粒子 点 间 实 际 距离 ，4d 表示 理论 距离 ， 在 本 文中 理论 值 等 


正 。 首 先 需要 计算 质点 间 约 束 变化 最 快 的 方向 ， 即 距离 约束 函 


数 的 梯度 YCCpi,p,) ， 这 样 就 可 以 确保 变化 的 方向 ， 然 后 通 


最 后 使 约束 中 的 两 个 粒子 恢复 
总 结 以 上 ， 可 以 得 到 


过 在 该 方向 上 矫正 Ap 变化 量 ; 
到 约束 范围 内 。 


Cp+APD = C(PJ+VC(D:Ap=0 3) 


为 了 能 更 好 地 控制 在 变化 方向 上 的 幅度 , 可 以 引入 一 个 标 


量 4 . 为 了 进一步 能 体现 重量 对 变化 的 影响 , 使 用 重量 信息 1 
的 倒数 WW 来 对 函数 优化 ， 从 而 Ap 就 可 以 表示 为 
C(p) 
Ap=-4:w———— VC(p) (4) 
1YCCD) 


当 使 用 距离 约束 时 ,约束 函数 的 基数 会 变 为 2， 则 对 应 的 
粒子 约束 方程 就 变化 为 


(5) 
ou Cpi, p,) 


DW 
IvC(pi,p)F 


此 方法 在 顺序 矫正 距离 约束 时 存在 一 个 严重 影响 模拟 效 
率 的 问题 。 如 图 1 所 示 ， 在 顺序 矫正 的 过 程 中 ， 首 先 矫正 的 距 
离 约 束 为 (a)， 即 a-b; 紧 接 着 又 要 对 〈b) 进行 距离 约束 矫 
正 ， 即 b-d， 这 就 使 得 这 次 矫正 会 对 上 一 步 矫正 中 b 的 位 置 产 
影响 ， 使 得 (a) 中 a-b 再 次 不 满足 其 距离 约束 ， 导 致 约束 
震荡 。 这 样 的 反复 修改 问题 使 布料 的 收敛 速度 大 大 降低 ， 而 最 
终 的 收敛 则 取决 于 模拟 踪 代 矫 正 的 次 数 。 这 一 现象 减缓 了 布 
料 收 敛 的 速度 ， 降 低 了 模拟 的 效率 。 


a C a 必 a [9 


VC(pi,p;) 


红 


d d d 
(a) (b) Ce 
图 1 震荡 问题 演示 医 


在 实际 模拟 中 ， 该 方法 虽然 能 提高 一 定 程度 的 模拟 速度 ， 
但 是 存在 的 约束 震荡 问题 会 使 整体 的 收敛 速度 依赖 于 和 迭代 矫正 
的 次 数 ， 从 而 使 得 模拟 存在 很 大 的 不 稳定 性 。 层 次 化 方法 恰好 
能 解决 这 个 问题 ， 并 提高 布料 动态 模拟 的 效率 。 


2 ”传统 层次 化 方法 


2.1 基于 位 置 的 层次 化 动态 模拟 

由 于 收敛 取决 于 欠 代 次 数 ， 这 导致 了 动画 的 帧 率 和 模拟 效 
率 都 不 能 达到 理想 的 状况 。 为 了 解决 收敛 及 效率 问题 ， 
Muller 国 提出 了 基于 位 置 的 层次 化 动态 模拟 方法 (HPBD) 。 该 方 
法 首先 通过 挑选 一 部 分 核心 粒子 进行 简化 网 格 ， 构 造 简化 层次 
约束 网 格 ， 接 着 利用 PBD 方法 对 原始 网 格 进行 模拟 ， 最 后 通过 
约束 网 格 中 各 简化 层 对 原始 层 进 行 二 次 约束 矫正 ， 在 整体 的 效 


间 
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果 上 得 到 了 进一步 的 提升 。 


在 使 用 基于 位 置 的 层次 化 模拟 方法 时 , 不 仅 需要 基本 的 位 


、 质 量 信息 1 以 及 各 种 约束 信息 ， 而 


且 还 需要 用 于 层次 网 格 控制 的 相关 变量 ， 如 层次 关联 信息 、 父 
言 息 些 


子 网 格 的 关联 信息 以 及 父子 网 格 中 粒子 间 的 对 应 关系 等 。 这 
言 息 如 表 1 所 示 。 
表 1 层次 化 额外 变量 
类 型 表示 符号 
粒子 类 型 Fine/ Coarse 
粒子 1 的 预测 位 置 图 2 约束 构建 
prediet 2.3 ”约束 网 格 构造 方法 
无 7711e 点 的 融入 点 层次 化 动态 模拟 方法 简化 部 分 的 主要 目的 在 于 构建 多 网 
格 约束 层次 系统 ， 使 得 系统 可 以 通过 逐 层次 、 不 同 程度 的 对 网 
粒子 2 ， J 间 的 距离 林 格 质 点 间 的 距离 约束 进行 矫正 。 这 样 一 来 ， 单 纯 地 依靠 原始 网 
" 格 而 产生 的 约束 震荡 现象 就 会 大 大 减 小 。 
保存 粒子 站 的 邻接 信息 i.Neighbors 每 个 层次 的 构建 包括 选择 原始 粒子 的 融入 粒子 、 各 层次 约 
粒子 二 的 父亲 节点 i.Parents 束 的 构建 以 及 确定 位 于 目标 层次 简化 网 格 的 Fine 点 。 具 体 构 
建 步骤 如 下 : 
人 Ws a) 要 根据 当前 Jevel 层 的 网 格 选择 在 Jevel 十 1 层 用 来 构 


> 层次 level 建 约束 网 格 的 粒子 ， 即 挑选 Jlevel 层 中 的 Coarse 点 ， 该 操作 
其 中 粒子 类 型 信息 表示 层次 化 过 程 中 粒子 所 处 的 层次 是 等 价 于 先 挑选 Fine 点 ， 然 后 剔除 。 其 流程 大 体 如 下 : 首先 判 
否 包 含 该 粒子 ， 即 当前 粒子 是 否 可 以 在 对 应 约束 层 中 剔除 ; 断 当 前 粒子 i.Neighbors 的 容量 是 否 大 于 大 〈 在 该 系统 中 天 
Fine 类 型 表示 在 目标 层 约束 网 格 中 将 不 再 包含 该 粒子 ， 如 图 ， 取 2)， 如 果 大 于 且 其 邻接 的 Fine 点 的 fine.Neighbors 也 
2 中 的 空心 点 ， 相 对 的 Coarse 类 型 则 是 目标 层 约束 网 格 中 包 ”大 于 到 -1 工 , 则 标记 为 可 剔除 的 Fine 点 ， 然 后 与 之 相关 联 的 粒 
子 7 从 其 的 j.Neighbors 剔除 粒子 ; 。 换 言 之 ， 就 是 如 果 粒 
括 的 粒子 ,如 图 2 中 的 实心 黑 点 ， 下 0e 点 的 融入 点 Pj 表示 目 。 地 j 与 有 关联 ， 则 从 其 邻 楼 信息 j.Neighbors 中 删除 粒子 
县 约束 网 格 中 将 该 Fine 点 的 信息 保存 的 点 ， 如 图 2 中 的 1 ”i。 最 后 剩 下 的 就 是 当前 层 挑选 出 来 的 Coarse 粒子 。 接 下 来 


多 | 


人 
en 


或 j2 点 ;dj 表示 粒子 i 与 j 间 的 距离 ， 即 两 粒子 间 的 弹 得 长 ”需要 继续 为 挑选 出 来 的 Fine 粒子 1 ( pi ) 添加 父亲 节点 


度 ; i.Neighbors 表示 当前 表示 粒子 i 的 临 接 粒子 ， 即 所 有 ”i.Parents .具体 流程 如 下 : 通过 对 i.Neighbors 进行 遍历 ， 
与 粒子 1 之 间 有 弹 得 连接 的 粒子 ，1i.Parents 表示 粒子 1 融入 如 果 对 应 的 粒子 类 型 是 Coarse , 则 将 其 填 入 到 i.Parents 即 


可 。 

Dj 后 的 Coarse 点 ， 如 图 2 中 的 fine 粒子 如 果 融 入 jl， 那 

么 其 父 节点 就 包含 j1 节点 , 在 这 个 局 部 它 融 入 了 j1, 下 个 局 部 

可 能 就 用 其 他 Fine 粒 子 也 融入 了 j1, 所 以 i.Parents 不 是 音 

一 元 素 ， 而 是 个 集合 ;Wi 表示 父亲 节点 (1.Parents ) 对 粒 

子 i 的 约束 矫正 影响 什 ， 计 算 方法 见 后 文 ，level 表示 在 层次 0 

约束 网 格 中 的 层次 位 置 ， 如 原 网 格 的 层次 数 为 6， 精简 n 次 则 b) 为 当前 level 层 的 生成 对 应 的 层次 约束 集合 ， 首 先 利用 

将 * 生 n 层 约 来 网 烙 ， 为 了 表示 当前 澡 作 的 层次 ,使 用 level 已 Me 类 型 粒子 ;找到 其 将 要 融入 的 粒子 j。 需 要 注意 的 是 ， 

表示 处 理 的 层次 。 j 的 候选 点 有 多 个 ,如 图 2 中 的 jl 与 j2， 如 果 jl 被 选 做 当前 
居 的 Coarse 点 时 ， Fine 类 型 粒子 与 Coarse 类 型 粒子 jl 


间 有 约束 ， 如 图 3 (a) 中 p; 与 Pj; 间 的 虚线 ， 删除 即 可 。 如 果 


201804.02137V1 


chinaXiv 


录用 稿 


Vine 类 型 粒子 邻接 的 粒子 与 融入 点 Pj; 无 连接 ， 此 时 需要 先 


删除 其 与 p; 的 连接 ， 然 后 添加 一 个 新 约束 ， 如 图 3 (b) 中 的 


DP; 。 如 果 Fine 类 型 粒子 邻接 的 粒子 与 融入 粒子 ;有 共同 的 


邻接 粒子 ， 删 除 此 约束 即 可 。 以 此 类 推 ， 剩 下 的 约束 即 当前 简 
化 层 的 约束 ， 至 此 当前 层次 的 简化 网 格 就 建立 完成 了 。 

c) 将 简化 后 的 约束 县 加 在 原 网 格 上 进行 收敛 加 速 。 这 部 分 
的 过 程 可 如 算法 2 所 示 。 

算法 2 HPBD 新 层次 构造 算法 
输入 : 上 一 层 网 格 约束 集合 。 
输出 : 新 层 网 格 约束 集合 。 

1. ”复制 所 有 层次 Jevel -1 的 约束 


2. for all 粒子 i do 

3. ifp,type= fine then 

4. for all p,.neighbors 1 do 
5. for all p ;neighbors m do 
6 if pineighbors,=p ;neighbors,, 
gp | 除 该 约束 

8. else 

9. 更 新 该 约束 

10. end if 

11. end for 

12. end for 

13. end if 

14. end for 


此 中 HPBD 约束 矫正 从 最 精简 层 开始 计算 ， 循 环 迄 代 至 最 
精密 层 ， 这 样 通过 父 节点 信息 与 子 节点 信息 的 传递 保存 了 布料 
的 关键 特征 ， 再 通过 利用 简化 网 格 约束 对 原 网 格 矫正 。 这 样 一 
来 ， 在 获得 信息 对 等 性 的 同时 不 仅 将 布料 信息 还 原 于 精密 原始 
网 格 ， 而 且 使 得 计算 量 大 幅度 减少 ， 保 证 了 效率 ， 并 且 保 证 了 
布料 模拟 的 真实 性 。 在 层次 信息 传递 的 步 又 上 ， 数 值 更 新 尤为 
重要 ， 具 体 计算 如 式 6) 所 示 。 


Xi Xit+ Dw (pj-q;) (6) 
eh 


其 中 :五 是 粒子 的 父亲 节点 的 索引 ， Wi 表示 忆 的 父亲 节 


点 ;对 其 的 影响 权重 ，p; 表示 粒子 i 的 父亲 节点 的 实际 位 


置 ，97 表示 粒子 i 的 父亲 节点 的 预测 位 置 。 关 于 Wy 的 计 


算 ， 本 文 使 用 粒子 i 与 当前 父 节点 p; 间 的 距离 di 与 粒子 i 与 


所 有 父 节 点 Pj; 间 的 总 距离 进行 求 比 进行 获取 ， 如 式 (7) 


《8) 所 示 。 先 求 取 未 单位 向 量化 的 值 W， ， 最 后 通过 向 量化 


处 理 得 到 权重 。 


CR 
” dj;/max(d;)+e 


去 a (8) 
Wi =W, /2 W, 


在 式 (7) 中 ，& 表示 一 个 特别 小 的 常数 ，max(d; ) 表示 


在 粒子 i 与 其 融入 粒子 j 之 间 最 大 的 距离 。 需 要 注意 的 是 ， 系 


统 中 dj 与 W; 是 随 所 在 层次 level 的 变化 而 变化 的 , 每 个 层次 


中 有 着 不 同 的 Coarse 粒子 ， 也 对 应 着 不 同 的 dj 与 W;。 


3 ”规则 网 格 层次 化 方法 


3.1 传统 方法 的 问题 
在 传统 方法 中 , 如 果 将 布料 网 格 不 同 层次 的 约束 网 格 提取 
出 来 进行 对 比 就 会 发 现 ， 随 着 层次 数 的 不 断 增加 ， 其 中 的 不 规 


则 多 边 形 数量 也 越 来 越 多 ， 如 图 4 所 示 。 


Bol 


(b) 第 三 


图 4 传统 层 


学 


从 以 上 可 以 看 出 ， 在 随 着 层次 数 的 增多 ， 传 统 算法 在 产生 


(c) 第 五 


MU 
S 
小 
Ct 
pl 


化 约束 网 格 


约束 网 格 的 过 程 中 存在 一 定 的 问题 ， 使 得 在 构造 层次 化 网 格 系 
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统 时 将 产生 大 面积 的 不 规则 多 边 形 ， 
影响 矫正 方向 ， 造 成 了 整体 的 约束 效果 与 


使 得 约 


束 网 格 的 质量 严重 


估 的 区 别 ， 这 将 导致 布料 模拟 的 最 终 | 


一 定 的 误差 。 


统 ， 但 是 会 存在 


接近 ， 但 在 某 些 特定 的 


区 域 中 ， 这 将 


然 文献 [11] 中 提 到 


传统 层次 化 方法 中 ， 在 约束 更 新 时 构造 属 了 
量 的 不 规则 多 边 形 ， 


里 论 状态 产 归 


性 敛 状 


hs | 
[一 


虽然 这 村 


的 随机 性 和 不 可 控 怕 


应 | 


法 中 ,层次 网 格 建立 世 


特别 多 的 资源 。 规 则 网 格 层次 化 方 济 


成 为 难题 ， 寺 
了 间 复杂 度 依赖 于 采 检 
和 这 两 种 算法 都 是 高 时 间 复杂 度 的， 所 以 在 建立 系统 上 会 花费 


使 得 方法 失去 可 靠 性 。 虽 


E 不 可 预 
与 原始 模拟 存在 


的 网 格 系 
在 整体 的 效果 


一 | 


的 采样 化 方式 可 以 解决 这 个 问题 ， 但 是 采样 


F 且 在 采样 层次 化 广 
和 网 格 重建 算法 ， 


将 利用 数学 构造 方 泣 


到 与 理论 效果 


得 到 了 质 的 飞跃 。 
3.2 ”层次 约束 网 格 的 建立 
对 于 层次 约束 系统 的 建立 


高 度 可 控 的 网 格 系统 ， 实 


约束 模型 ; 最 后 根 
在 挑选 次 层 
的 列 和 行 都 需要 不 同 的 


数 ， 例 如 对 网 格 


| 
4 


先 将 原始 网 格 作为 基准 
分 析 ， 挑 选 出 最 适合 的 决策 函数 ， 从 而 确定 模拟 材料 的 的 
8 的 Fine 点 进行 构造 层次 约束 网 格 。 
粒子 (Coarse 点 ) 的 过 程 需 要 一 定 的 坟 
函数 f(X) 和 8&(y) 对 


各 自 点 的 个 数 进行 确定 ， 然 后 再 通过 离 


首先 需要 确定 者 


; 然后 通过 对 者 


岗 局 部 与 整体 的 双重 规则 的 要 求 ， 达 
且 在 构建 层次 约束 网 格 的 效率 上 


实现 


入 
漆 
路 
这 
王 
奖 
< 


郑 华 ， 等 : 


再 比如 , 当 行 的 粒子 初始 个 数 为 100, 列 的 初始 粒子 为 30， 
希望 在 此 网 格 中 比较 密集 的 粒子 能 够 相对 绥 慢 的 递减 ， 比 较 疏 
松 的 列 粒子 能 够 均匀 递减 ， 所 以 此 时 在 列 的 简化 过 程 中 ， 本 文 
使 用 线性 函数 作为 决策 函数 ， 在 行 的 简化 过 程 中 ， 使 用 余弦 函 
数 作为 决策 函数 。 


f (1)=Ek cos(@l) +b, 


(11) 


g(1,)=kl+b, - 


在 式 (11) (12) 中 : 有 、k、0@9、b 和 bb, 均 为 常数 ; 


了 (4,) 为 规则 网 格 的 行 粒子 个 数 ，8 (4 ) 为 规则 网 格 的 列 粒子 


数目 ，! 代表 处 理 的 所 处 层次 。 


准 网 格 。 本文 首 


准 网 格 的 特点 进行 


层次 


R 策 函 


的 具体 位 置 。 如 
过 程 中 ， 列 和 行 都 使 用 
是 属于 递减 类 型 区 


网 格 往往 需要 根 


例如 在 图 5 中 ， 


粒子 个 数 也 为 5， 


次 约束 网 格 ， 使 得 其 变化 保持 一 致 ， 所 以 本 文 使 用 
函数 进行 构造 约束 网 格 ， 保 订 


(a) 原始 


Cl 


的 行 粒子 个 数 ， 8 (4 ) 为 规则 网 格 的 列 粒子 数目 ; 


的 所 处 层次 。 


f(1)=kI+Dh 


g(1,)=kI+h, 


艇 的 方式 确定 每 个 粒子 
if 示 ， 在 图 5(a) 到 (b), 最 后 简化 到 (d) 的 
了 相同 的 线性 函数 。 在 选择 函数 上 ， 凡 
数 都 可 以 作为 约束 网 格 的 区 
时 实际 情况 确定 其 在 行列 上 不 同 的 决策 函数 。 
首先 在 网 格 行 的 初始 粒子 个 数 为 5， 列 的 初始 
且 在 常规 的 减 函数 中 希望 存在 比较 多 的 


策 函 数 。 不 同 的 


Ml 


相同 的 线性 
FE 行 与 列 的 变化 一 致 性 。 
(9) 


(10) 


(c) 第 二 层 


灵 次 化 网 格 构造 


(d) 第 三 


总 


在 式 (9) (10) 中 : 后 、b 均 为 常数 ， 了 (7 ) 为 规则 网 格 


! 代表 处 理 


在 这 个 过 程 中 , 首先 需要 确定 最 小 层次 中 行列 的 粒子 个 数 
和 最 大 层次 中 行列 的 粒子 个 数 ， 如 图 5 中 最 大 层次 的 行列 粒子 
个 数 均 为 2， 最 小 层次 的 行列 个 数 均 为 5。 在 此 例 中 ， 最 小 层次 
的 行 粒子 个 数 为 100, 最 大 层 的 行 粒子 个 数 为 2。 最 小 层次 的 列 
粒子 个 数 为 30， 最 大 层 的 行 粒子 个 数 为 2。 在 此 便 可 以 确定 其 
中 变量 的 具体 值 ， 从 而 可 以 确定 出 各 层次 的 粒子 具体 位 置 。 本 
例如 表 2 所 示 。 


表 2 粒子 各 层次 变化 


变化 ”初始 第 1 第 2 名 3 名 4 ”中 间 最 后 

刁 层 属 导 E: 层 
行 100 85 72 53 34 人 2 
列 30 26 22 8 4 2 2 


在 确定 每 层 中 每 行 的 粒子 个 数 与 每 列 的 粒子 个 数 后 , 需要 
通过 离散 的 方式 构造 网 格 中 每 个 点 的 位 置信 息 。 例如 图 5(a) 中 
的 粒子 通过 g(y) 确定 了 下 层 每 行 的 粒子 个 数 ， 同 样 通过 
了 (x) 确定 了 每 列 的 粒子 个 数 ， 然 后 计算 出 最 长 行 ( 列 ) 端点 
之 间 的 距离 ， 求 其 离散 增长 步 长 ， 从 而 确定 出 下 一 层次 的 网 格 
模型 (图 5 (b) )。 然 后 图 5 (b) 重复 同样 的 操作 ， 直 至 到 达 定 义 
的 终止 条 件 。 在 这 里 ， 称 图 5 (b) 为 (a) 的 父 网 格 ， 称 图 5 (c) 
为 (b) 的 父 网 格 ， 并 以 此 类 推 。 


此 时 , 父 网 格 的 运动 轨迹 如 何 确定 又 成 了 另 一 个 核心 问题 。 
笔者 使 用 了 一 个 传统 的 解决 方案 ,利用 最 短 距离 确定 其 同位 素 ， 
即 在 构建 层次 系统 时 ， 距 离 该 粒子 最 近 的 原 网 格 粒子 将 与 该 粒 
子 有 相同 的 位 置 与 位 移 关 系 ， 这 就 确保 了 约束 网 格 的 运动 特征 
与 原 网 格 的 一 致 性 。 


3.3 ”模拟 关联 

确定 了 层次 系统 ， 如 何 将 约束 信息 串联 起 来 是 这 部 分 的 核 
心 问 题 。 正 如 图 6 所 示 ， 该 部 分 为 图 5(b) 到 (c)， 其 中 蓝 色 的 
为 上 层 删 选 掉 的 粒子 ， 即 Fine 点 ， 黑 色 的 为 本 层 的 粒子 ， 即 
Coarse 点 ( 见 电 子 版 ), 并 且 也 可 以 将 (c) 视 为 (b) 的 父 网 格 。 
为 了 将 父 网 格 的 约束 矫正 信息 能 够 传递 给 子 网 格 ， 需 要 在 构建 
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录用 稿 


约束 网 格 系统 时 ， 确 定子 网 格 粒子 在 父 网 格 的 父亲 粒子 们 ， 并 
确定 对 应 的 权 值 信息 。 


[ 亚 


(b) 父子 粒子 关联 


图 6 父子 粒子 关系 
对 于 图 6 (b) 中 的 4、Bb、C 与 P 均 存在 于 三 维 空间 ， 则 
必然 存在 一 组 一 维 变量 W 、W2 、Ws， 使 得 
(13) 
P=wat+wbt+wc 
这 个 式 子 等 价 于 
Px Ul 已 。 Clx Wi 
Py = 0 b, C1y Ww 
Pp: di, 已 Cl1z Ws 
假如 使 用 P 表示 Fine 粒子 的 位 置 矩 阵 ，A4 表示 其 父亲 
节点 的 位 置 乍 阵 ， 人 三 表示 其 权重 矩阵 ， 则 有 
P=A-:W (14) 
则 可 以 推导 出 
W=A1.P (15) 


求 出 父子 权 值 后 ， 便 可 以 在 模拟 时 使 用 对 应 权 值 更 新 子 节 


点 位 置 ， 从 而 避免 
3.4 ”模拟 流程 


了 约束 震 落 问 题 ， 加 快 收敛 速度 


对 于 本 文 的 
层次 系统 ， 其 次 在 


层次 化 方法 ， 在 模拟 的 过 程 中 


o 


首先 需要 确定 


约束 矫正 的 过 程 中 


初始 
的 位 置 ， 然 后 在 约 
图 5 中 利用 图 
格 (a) ,这 样 的 层次 
矫正 中 使 用 的 迭代 
中 ， 子 网 格 中 的 粒 


的 信息 ， 正 如 第 3 章 中 的 式 (6) 所 示 。 整 体 流程 如 下 算法 3 


不 。 


算法 3 本 文 
布料 网 格 ， 
布料 动态 运 


输入 : 
输出 : 


5 (qd) 矫 正 (c)， 利用 
化 矫正 可 以 减少 原始 动态 模拟 方法 中 在 约束 
在 具体 的 矫正 过 程 
父 网 格 中 父亲 节点 们 获取 其 新 的 位 置 


首先 需要 对 原始 布料 进行 


化 模拟 ， 即 通过 力学 方程 预测 出 下 一 时 间 步 长 布料 所 有 点 


束 网 格 中 从 最 粗 的 一 层 开始 回归 


矫 


正 ， 即 如 


至 


(c) 矫正 (b)， 直 


次 数 和 总 体 的 约束 数量 。 
子 利用 


层次 化 模拟 管线 
外 界 环境 参数 (如 重力 、 
动 轨迹 。 


1. ”构造 约束 网 格 
for all 粒子 i do 


初始 化 和 


end for 


=x,V, =V ,Ww =1/m, 


i 


for all 粒子 i do 


2 
3 
4 
5 loop 
6 
2 
8 end for 


Vv, V+At:w,:f,,() 


到 原始 网 


所 
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9 for all 粒子 i do 
10 Pp; xX,+At:v, 
11 end for 
12 loop 和 迭 代 次 数 
13 原 约束 矫正 
14 end loop 
1 for all 约束 网 格 (逆序 ) do 
16 Pp; -Pit 六 W(X 一 Xi predict ) 
jei.Parents 
17 回归 约束 矫正 
18 end for 
19 for all 粒子 i do 
20 vv €(p,—xX)/At 
2 AD 
22 end for 
23 end loop 
在 此 处 ， 为 了 防止 预测 位 置 表示 符号 的 重复 带 来 的 混淆 ， 
其 本 的 位 置 预测 使 用 p; 表示 , 父 类 的 位 置 预 测 使 用 XX; wredier 表 
示 。 与 原 层次 方法 不 同 的 是 ， 在 回归 约束 矫正 的 过 程 中 ， 不 再 
是 将 父 层 的 矫正 信息 传递 至 子 层 ， 而 是 直接 传递 到 原始 层 ， 这 
样 不 再 过 分 依赖 某 个 关键 粒子 ， 便 使 得 整体 的 效果 接近 。 
4 ”实验 结果 分 析 与 效果 对 比 
4.1 添加 风 场 
为 了 校 验 基 于 规则 网 格 的 层次 化 动态 模拟 方法 ， 本 文采 用 
文献 [13,14] 提 出 的 随机 风 场 模型 ， 并 结合 Perlin 噪声 函数 产 
生 的 钟 摆 05 风 场 进 行 测试 。 首 先 使 用 三 角形 网 格 来 表示 风 场 中 
布料 模型 ,在 分 析 过 布料 的 网 格 结构 后 , 可 以 确定 每 个 三 角形 的 
顶点 也 是 其 相 邻 三 角形 的 顶点 。 在 计算 风力 时 ， 由 于 使 用 的 模 
型 是 离散 型 ,所 以 先 计算 顶点 在 每 个 三 角形 面 片上 的 风力 , 再 通 
过 加 权 平 均 来 求 取 每 个 顶点 所 受 的 风力 ， 最 后 根据 空气 动力 学 
里 论 中 空气 流体 对 物体 的 作用 力 来 计算 具体 三 角形 面 片 所 受 的 
风力 。 计 算 方法 如 式 〈16) 所 示 。 
(16) 


1 2 2 
F = 一 CDV 9， 
i 2 Pp i 


其 中 : CC 表示 风力 系数 ， 一 般 取 值 在 0.3 一 0. 6 间 ; P 表示 空 


气 密度 ， 大 小 一 般 取 2.37e 一 3g /11 ; Y 表示 风速 ，S, 表 示 


力 。 对 了 


角形 盏 


片 的 


片 的 受 力 ， 再 通过 求 


4 


于 采样 的 
行 模拟 。 对 了 


和 


.2 ”模拟 效果 对 比 


离散 点 上 的 力 ， 需 要 
取 平 均值 来 计算 。 


本 文 使 用 ] 


层次 化 模拟 方法 与 规则 网 格 
本 文 方法 ， 首 9 


积 ; 而 五 表 示 第 i 个 三 


形 面 片 所 受 的 
E 求 取 所 有 与 其 相关 联 的 三 角 面 


基于 位 置 的 动态 模拟 、 传 统 的 层次 化 方法 、 基 


上 用 自 定 义 校 准 


函数 确定 


层次 化 方法 对 布料 
E 将 原始 布料 视 为 基准 网 格 ， 然 后 


网 格 进 


层次 化 约束 系统 ， 紧 接着 利用 


层次 间 
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录用 稿 郑 华 ， 等 ; 基于 规则 网 格 和 料 动态 模拟 


关系 确定 粒子 之 间 的 关联 ， 使 得 约束 信息 的 有 效 传递 ， 最 终 完 度 的 约束 矫正 ， 使 得 模拟 效果 远 好 于 仅仅 增加 和 迭代 次 数 、 增 加 


整 建立 层次 化 系统 。 在 模拟 的 过 程 中 ， 首 先 需要 对 原始 网 格 利 约束 矫正 次 数 的 原 模拟 方法 。 实 验 模拟 效果 如 图 7 所 示 。 
用 动力 学 方程 进行 初步 的 位 置 预 估 ， 然 后 从 最 精简 层 到 最 精 台 在 以 下 的 测试 中 ， 核 心 硬件 架构 为 x86， 处 理 器 为 Intel 
层 分 别 利用 层次 系统 中 对 应 的 层次 约束 进行 矫正 ， 使 得 布料 网 。 酷 害 (型 号 i7-4710MQ) ,搭载 Windows 系统 ， 同 时 使 用 开源 图 
格 信息 在 层次 传递 下 保持 一 致 。 形 库 0penGL 使 用 单线 程 绘 图 方式 进行 泻 染 。 

在 基于 位 置 的 动态 模拟 中 ， 为 了 能 更 好 地 达到 收敛 效果 ， 首先 对 基于 位 置 的 动态 模拟 方法 进行 模拟 采样 , 然后 在 同 
只 能 单一 地 增加 约束 的 矫正 次 数 ,存在 着 大 量 的 约束 震荡 问题 。 ” 样 的 情况 下 分 别 对 传统 层次 化 方法 、 基 于 采样 的 层次 化 模拟 方 


在 层次 化 的 动态 模拟 中 ， 虽 然 添 加 了 更 多 的 约束 , 但 是 由 ”法 、 基 于 规则 网 格 的 层次 化 方法 进行 模拟 采样 。 以 下 是 在 重力 
于 不 同 层次 约束 的 添加 却 使 得 迭代 次 数 的 减少 ， 而 且 在 总 的 约 ”， 为 0. 98 牛顿 、 迁 代 次 数 为 5、 时 间 为 3 s 的 标准 条 件 下 进行 自 
束 数量 上 比 原 模拟 少 了 许多 ， 这 就 直接 造成 了 效率 的 提升 。 落体 实验 的 对 比 图 (图 7- 第 一 行 )。 
且 在 模拟 的 效果 上 ， 由 于 不 同 层 分 别 对 整体 各 范围 进行 不 同 程 


we 


原始 模拟 传统 层次 化 方法 采样 层次 化 方法 本 文 层次 化 方法 


原 方法 (a) 原 方法 (b) 原 方法 (c) 原 方法 (d) 
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为 了 测试 布料 在 风 场 中 形变 的 稳定 程度 , 本 文 在 定向 钟 摆 


风 场 下 进行 测试 ， 效 果 如 民 
力 为 0.98 牛顿 ， 和 迭代 次 数 为 5， 约束 网 格 
完成 。 其 中 上 方 一 行为 本 文 方法 下 的 模拟 效果 


7 第 2、3 行 所 示 。 该 部 分 实验 在 重 
层次 数 为 3 的 情况 下 


统 方法 的 模拟 效果 。 
收敛 速度 ， 但 是 本 文 方法 石 


通过 对 比 可 以 看 出 ， 传 统 方法 虽然 提高 了 
E 动 态 条 件 下 正如 在 风 场 下 的 效果 所 


F 面 一 行为 传 


示 ， 原 始 方法 下 布料 的 拉 伸 状态 与 
与 


明基 于 规则 网 格 的 


原 网 格 差距 较 大 。 这 巧合 说 


等 方面 上 做 得 更 加 突 时 
为 了 将 各 种 方法 


理论 状态 更 加 接近 ， 在 稳定 性 


的 性 能 以 及 | 


H， 在 粒子 位 置 的 收敛 速度 上 效果 明显 。 
改 敛 速度 显 像 地 表示 出 来 , 笔 


者 使 用 最 常用 的 弹 筑 变化 率 对 对 传统 方法 、 原 始 基于 位 置 的 动 


态 方法 (PBD)、 原 始 


以 及 基于 规则 网 格 


层次 化 方法 (HPBD )、 采 样 层 次 化 方法 (SHPBD) 
的 层次 化 方法 (GHPBD) 进行 测试 比 对 ， 以 说 


明 各 种 方法 最 终 收敛 状态 之 间 的 差距 。 测 试 数据 如 表 3 所 示 。 
表 3 稳定 弹簧 相对 变化 率 比 较 


方法 HPBD SHPBD GHPBD 
层次 化 1 层 63% 46% 37% 
层次 化 3 层 48% 34% 23% 
层次 化 5 层 41% 26% 19% 
而 在 传统 方法 中 ， 最 终 收 化 状态 下 的 相对 变化 率 为 52%， 


的 时 间 为 42 ms。 


虽然 在 相对 效果 好 于 最 原始 的 
的 耗 时 却 远大 于 PBD 系列 方法 。 


PBD 方法 ， 


但 是 在 模拟 的 过 程 中 


例如 在 单 帧 模拟 耗 时 上 ， 传 统 


方法 耗 时 为 132 ms 。 而 淄 


基于 采样 的 方法 分 别 为 97 ms、34 ms 利 


通过 在 相同 条 件 下 对 原始 动态 模拟 方法 、 传 统 的 层次 化 动 


态 模 拟 方法 与 基于 规则 网 格 包 


收敛 ， 并 且 在 布料 的 动态 
4.3 ”收敛 速度 对 比 

布料 网 格 随 着 时 间 
能 够 说 明 一 种 方法 的 


本 文 提出 的 基于 规 贝 


沪 


的 动态 模拟 、 传 统 层次 方法 和 
0 45 ms， 而 本 文 方法 


郑 ” 华 ， 等 :基于 规则 内 格 的 层次 化 机 入 动态 模拟 方法 


图 7 模拟 效果 对 比 


一 一 本 文 方法 。……. 传统 方法 “一 一 一 原 哲 模 拟 


8 ”弹簧 变化 率 对 比 


通过 图 8 可 以 看 出 ， 使 用 本 文 提 出 的 基于 位 置 的 层次 化 模 


拟 方法 时 ， 虽 然 耗 时 相 比较 略微 增长 ， 但 是 布料 的 整体 弹簧 变 
化 率 能 更 快 地 趋向 于 布料 整体 收敛 位 置 。 


结束 语 


本 文 首先 对 布料 分 别 利用 基于 位 置 的 动态 模拟 方法 、 传 统 


层次 化 方法 模拟 方法 和 本 文 提出 的 基于 规则 网 格 的 层次 化 方法 
进行 模拟 ， 然 后 通过 测试 说 明了 基于 规则 网 格 层 次 化 方法 在 动 
条 件 下 的 稳定 性 ， 最 后 通过 对 比 说 明了 本 文 方法 在 效率 及 效 


一 


面 的 优越 性 。 在 传统 方法 的 实现 过 程 中 ， 随 着 层次 数 的 增 


加 速 约束 网 格 中 存在 的 不 规则 多 边 形 越 来 越 多 ， 这 一 网 格 


I 质量 问题 使 得 布料 的 收敛 趋 于 未 知 ， 导 致 了 部 分 区 域 的 收敛 


里 论 上 存在 一 定 的 差异 ， 


且 在 构造 层次 系统 的 过 程 中 对 于 


层次 化 模拟 方 


法 的 模拟 对 比 可 知 


化 方法 ， 布 料 的 收敛 状态 更 快 


地 接近 理想 状态 ， 证 


车 


模拟 


的 层次 化 模拟 方法 确实 加 速 J 
效果 稳定 。 


j 位 置 趋向 于 最 终 位 置 的 变化 ， 


传统 层次 化 方法 和 本 文 方法 使 ) 
并 从 这 个 过 程 中 比较 


在 此 ， 笔 者 分 别 对 原始 方法 、 
同样 的 测试 条 件 进行 逐 帧 比 对 ， 
广 敛 速度 。 


通过 对 比 发 现 随 着 时 间 的 移动 ， 层 次 化 模拟 的 优势 逐渐 显 


现 出 来 ， 传 统 方法 在 中 


本 现 出 优势 ， 但 是 后 期 与 基于 位 置 的 


动态 方法 表现 收敛 速度 开始 


显示 出 其 优越 性 ， 


果 持 领 


法 的 收敛 加 速效 果 , 测 试 使 月 
时 的 布料 形变 指数 ) 逐 帧 比 对 。 


下 降 ， 而 本 文 方法 从 中 前 期 就 开始 
E。 为 了 更 好 地 说 明 本 文 方 
日 了 弹簧 变化 率 (表示 布料 趋 于 稳定 


人 No 


而 对 于 采样 方法 , 由 于 依赖 于 采 
样 算法 与 重 构 算法 ,使 得 其 该 部 分 的 测试 数值 局 限于 一 个 范围 ， 
所 以 本 文 不 提供 其 该 部 分 的 弹 


率 测试 数据 。 测 试 如 图 (8) 所 


季 的 控制 率 比 较 低 。 传 统 层 次 化 方法 为 了 解决 最 初 的 约束 震 
荡 问 题 而 带 来 的 细节 控制 问题 一 直 未 被 解决 ， 本 文 针 对 这 个 问 
题 提 出 了 基于 规则 网 格 层次 化 方法 ， 创 新 性 地 使 用 行列 决策 函 
数 来 确定 层次 约束 网 格 ， 并 结合 采样 层次 化 方法 中 的 关联 方式 
成 功 解决 了 传统 方法 的 效率 问题 、 约 束 震 荡 问 题 以 及 约束 网 格 
量 问题 导致 的 收敛 不 确定 问题 , 在 父子 粒子 的 关系 控制 上 达 
到 新 的 高 度 ， 在 模拟 细节 上 有 了 更 好 的 表现 ， 在 风 场 的 动态 测 
试 中 表现 良好 。 
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